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>< ■ Abstract 

We calculate the inelastic scattering cross sections to populate one- and 
two-phonon states in heavy ion collisions with both Coulomb and nuclear ex- 
citations. Starting from a microscopic approach based on RPA, we go beyond 
it in order to treat anharmonicities and non-linear terms in the exciting field. 
These anharmonicities and non-linearities are shown to have important ef- 
fects on the cross sections both in the low energy part of the spectrum and 
in the energy region of the Double Giant Quadrupole Resonance. By prop- 
erly introducing an optical potential the inelastic cross section is calculated 
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semiclassically by integrating the excitation probability over all impact pa- 
rameters. A satisfactory agreement with the experimental results is obtained. 
PACS : 21.60. Ev; 21.60.Jz; 24.30. Cz; 25.55 Ci; 25.70.De 

Keywords: Coulomb and nuclear excitation, multiphonon states, anharmonicity and 
non-linearity in RPA, heavy ion collisions. 

I. INTRODUCTION 

All theoretical approaches used to calculate the cross section for the multiple excitation 
of Giant Resonances (GR) in heavy ion collisions are based on a semiclassical description 
of the process [|l|, where the excitation of one reaction partner is assumed to be due to the 
action of the mean field of the other and is treated quantum mechanically while the relative 
motion is determined classically. 

For each eigenstate a of the internal hamiltonian of one nucleus, one can calculate its 
excitation probability by perturbation theory or by solving a system of coupled equa- 

tions. This is done by integrating the equations of motion along the classical relative motion 
trajectory corresponding to the impact parameter h. The total excitation cross section cTq, is 
then evaluated by integrating the probability over all the impact parameters starting from a 
minimum one hmin- In Coulomb excitation studies the value of the latter is chosen according 
to a systematics following some prescription based on the condition that the contribu- 
tions from the nuclear field should be eliminated. Even so, however, some ambiguities are 
present since the calculated cross sections can vary appreciably for small variations of hmin- 
Moreover, when the bombarding energy is not very high and the two nuclei are not very 
heavy, the nuclear excitation is the dominant process. In this situation one cannot apply 
that procedure because in principle one should add more internal trajectory for the determi- 
nation of cTq. On the other hand, the trajectories corresponding to small impact parameters 
would not contribute too much to the inelastic cross section if the absorption due to all 
other channels is taken into account. This can be done by introducing an optical potential 
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as was already done in a qualitative way in ref. [Q. 

In this paper we present calculations for the excitation cross section of one- and two- 
phonon states in the ^'^Ca + ^°Ca reaction at 50 MeV/u for which experimental results 
exist [0. The calculations are done within the extended RPA model described in our previ- 
ous works where we have introduced anharmonicities in the internal hamiltonian and 
non-linear terms in the external field. This model has been successful in the description of 
the excitation of the double giant resonances, reducing the discrepancy between the mea- 
sured cross section and the standard theoretical estimate. Here the model is extended by 
introducing an optical potential in order to avoid the uncertainty on the integration over the 
impact parameter. Since the optical potential takes into account the absorption due to all 
channels, we have introduced a procedure in order to avoid to double count the effects of the 
channels explicitly included in our calculations. In the next section we will recall briefly our 
model and extensively describe its improvements. In section III we present our results and 
the quantitative comparison with the experimental findings. We then draw our conclusions 
and discuss some perspectives. 

II. THE MODEL 

The best microscopic theory to describe collective excitations in nuclei is the RPA whose 
hamiltonian can be written as 

Hrpa = Y.EuQIQu (2.1) 

V 

where the phonon creation operator is 

qI = Y.^x;^^1^-y;^b,^). (2.2) 

The bosonic operators B are the lowest order terms of the bosonic expansion of the fermionic 
operators 

aja, ^ + (1 - V2) 5: + ... (2.3) 
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Here, the index p (h) labels the particle (hole) states with respect to the Hartree-Fock ground 
state. The other terms after the first one correct for the Pauli principle. 

In the harmonic RPA hamiltonian only the Vph^p'h' and Vpp'^hh' terms of the residual 
interaction are taken into account. If we consider also the other terms lpp',p"p"'' Vhh',h"h"', 
Vpp',p"h and Vph^h',h" and introducing the mappings 0] 

alttp' — > {alap')B = J2h Bl^Bpih 

«/i4' — ' («/i4')b = EpBlf^Bph' 
one ends up with a hamiltonian containing cubic, quartic, etc, terms in the phonon creation 
and annihilation operators. In the space spanned by one- and two-phonon states the bosonic 
hamiltonian is 

h = y: eMQu + [ E yuLu^QlQlQr^. + E yuL.u,.,QlQlQM + h.c. (2.5) 

where V"^^ (V^^^) are the matrix elements connecting one- with two-phonon states (two- with 
two-phonon states). The eigenstates of the hamiltonian ( p.5|) are 

l*a >= Ec"k > + E €,^1^11^2 > (2.6) 

U V\V2 

and the corresponding eigenvalues do not form a harmonic spectrum. 

In the semiclassical models of grazing ion-ion collisions the excitation of one of the two 
nuclei is due to the mean field of the other. Since the mean field is a one-body operator, the 
excitation operator has the following form 

Wit) = E < a|f/B(R(t))|/? > 4^/3 (2.7) 

where Vb is the mean field of the other nucleus. The time dependence comes in through 
the relative distance R between the two nuclei. In the standard approach Wif) is linear in 
the phonon operators because only the ph terms of eq. (p.7| ) are considered and the lowest 
order boson expansion is taken. If we include also the pp and hh terms, their mapping (eq. 



2.4) leads to a quadratic form in Q 
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W = W'' + Y: Wl'Ql + h.c. + Y: Wll,QlQ,. + E WZQlQl + h.c. (2.8) 



The first term in eq. (|2.8|) represents the interaction of the two colliding nuclei in their ground 
state, in the present case it has also an imaginary part which describes the absorption due 
to the non elastic channels. The W^^ part connects states differing by one phonon, the W^^ 
term couples excited states with the same number of phonons, while allows transitions 
from the ground state to two-phonon states. All the form factors W are calculated by 
double-folding the Coulomb and nuclear nucleon-nucleon interactions with the Hartree-Fock 
ground state density of the projectile and with the ground state density or the transition 
densities of the considered excited states of the target. 

In the space of the ground state and the > states we can cast the Schrodinger 
equation into a set of linear differential coupled equations for the time dependent amplitude 
probabilities Aa{t). Then the cross section is calculated non-perturbatively as described 
in ref. where we integrated the probability of exciting the state |$q > starting from a 
minimum impact parameter. In the calculation presented here we integrate over all impact 
parameters since we have introduced in W^*^ the optical potential which, in an effective way, 
takes care of the most inner trajectories. 

The imaginary part Wim of the optical potential is usually determined by fitting the 
experimental elastic cross section. This potential describes the absorption due to all non- 
elastic channels. Therefore, it cannot be inserted directly in W^'^ (eq. (|2.7| )) since the 
absorption due to the inelastic channels explicitly included in the coupled equations would 
be counted twice. 

Let us first discuss how to solve this problem when no anharmonicities are present and 
therefore, the states |<I>q, > are pure multiphonon states. In such a case one can solve the 
Schrodinger equation in a semiclassical approach by integrating it along each classical relative 
motion trajectory. The state of the system |\E' > is a coherent state and the probability to 
excite rriv times a phonon is 



where is the average number of z^-phonons in >. In the above equation, as well as 
in the following discussion, the dependence on the impact parameter h is understood. The 
superscript "0" refers to the fact that only the absorption due to the multiple excitation of 
phonons is taken into account. In such a case 

Pi. = (2.10) 

where = Ny. We stress that the same survival probability of the ground state appears 
factor in all the probabihties in eq. ( p^ . 
The survival probability associated with the imaginary optical potential Wim is calculated 

as 

2 r+oo 

PZ. = exp W,m{t) dt} (2.11) 

where the integral is again done along a classical trajectory. The de-population of the ground 



state due only to the neglected channels can be in principle calculated as in eq. ( |2.11| ) but 
with an auxiliary imaginary potential W which does not contain the absorption due to the 
adopted ones. Then 

PZ. = Ks. X pZ. (2-12) 

and 

Pu,m,, = Pu,m^ ^ PZ- ■ (2.13) 

When anharmonicities are taken into account, the state of the system is no more a 
coherent state. The probability to excite the state |$a > is equal to 

po = \A^? (2.14) 

where is solution of the coupled equations of motion without any imaginary potential. 
Therefore, contains only the absorption due to all the adopted channels. The remaining 
part, due to all the other channels, can be introduced by writing, in analogy with eq. ( p. 131) , 
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P, = P°xF;,.. (2.15) 

This equation can be formally derived by assuming that the absorption due to the excluded 
channels is the same in all the adopted ones. This is certainly an approximation, however 
we would like to emphasize that many important inelastic channels are explicitly taken 
into account in the coupled equations and that we solve the latter exactly. Therefore, the 
corresponding absorption is calculated correctly, including the Q- value effects. The unknown 
auxiliary imaginary potential W can be eliminated by inserting eq. (|2.12|) in eq. (|2.15| ) 



pw 

P^ = P>^ (2-16) 
g-s. 

which is the expression we have used in order to calculate the inelastic cross section. We 
would like to stress that the part of the nuclear absorption that corresponds to non inelastic 
channels is often taken into account as a sharp cut off transmission coefficient. So the 
introduction of the imaginary potential can be seen as an important improvement. 



III. RESULTS AND DISCUSSION 

The above described model has been applied to the reaction ^''Ca on ^''Ca at E/u = 50 
MeV. The one-phonon basis has been obtained with a self-consistent HF+RPA calculation 
with Skyrme interaction SGII [l^. Only the most collective one-phonon states, exhausting 
at least 5% of the relevant EWSR, are taken into account. They are listed in table |. We then 
have considered all possible two-phonon states that can be constructed out of them, with 
all possible values of the total angular momentum L, and in this space we have diagonalised 
the hamiltonian ( |2.5| ) to get the states (p.6|) . In table H we have reported some properties 
of the quadrupole states, each one labeled with the name of its main component and whose 
unperturbed energy is given in the second column. In the third column there are the energy 
shifts due to the anharmonicities. Their overlaps with the single and double ISGQR states 
are shown in the last two columns. Similar tables for the GDR states are reported in 0. 
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The elementary nuclear form factors W to pure one- and two-phonon configurations 



(eq.(p.8[)) were calculated by double folding the M3Y nucleon-nucleon interaction with 
the RPA transition densities. The transition matrix elements between mixed states |$q, > 
were computed by mixing the elementary form factors according to the unitary transfor- 
mation ( |2.(j| ). The same procedure was used with the Coulomb interaction to calculate the 
Coulomb form factors. The relative motion trajectories were determined by solving the clas- 
sical equation of motion in the presence of both Coulomb field and real part of the nuclear 
potential. 

The real part of the optical potential was obtained by double folding the M3Y nucleon- 
nucleon potential with the Hartree-Fock densities of the two nuclei while its imaginary part 
was chosen with the same geometry and multiplied by a scale factor whose value (0.627) 
was determined by a fit to the experimental elastic cross section for the collision ^°Ca on 
^oCa at E/u = 50 MeV of ref. . 



In these calculations both the nuclear and Coulomb excitations were included. Actually, 
the Coulomb excitation alone does not produce a sizable cross section because the colliding 
nuclei are not very heavy, but when it is considered together with the nuclear excitation 
it produces an interference effect which can be important. This is due to the fact that on 
one hand we have a coupled channel effect and, on the other hand, some two-phonon states 
are excited only when both fields are acting. This was clearly demonstrated in our previous 
work [0. 

Since our calculations are based on a discrete RPA we get a discrete excitation spectrum 
and a cross section corresponding to each state >■ The energy differential cross 
sections presented in fig. |l| are obtained by summing up all the contributions coming from the 
states |$Q, > after a smoothing of each individual line by a Lorentzian with a 3 MeV width. 
The dashed line refers to a calculation where the internal hamiltonian is harmonic and the 
external field is linear. The solid line corresponds to a calculation where the anharmonicity 
and non-linearity were introduced, which produce a sizable increase with respect to the 
standard case. In the figure we can clearly distinguish three energy regions. The cross 
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sections given in tables |TT| to are obtained by summing up the a^s for the discrete states 



> lying in each region. As already observed in ref. I^Q? the increase at low energies 
is due both to the anharmonicities and non-linearities. In particular, the anharmonicities 
are important because the low lying two-phonon states can be excited by the W^^ part of 
the external field through their large one-phonon component. At high energies the main 
contribution comes from the non-linearities because their presence increases the number of 
excitation routes. This is seen better in table where the excitation cross section in the 
double giant quadrupole resonance energy region is reported. For each multipolarity we have 
summed the excitation cross section in the energy region between 28 and 38 MeV, and this is 
done for four different cases as shown in the table. The L=3 contribution is due to the HEOR 
at 31.33 MeV, while the L=0,2 and 4 contributions are dominated by the double excitation of 
the double ISGQR. As we can see in table 0, the non-linear terms are also responsible for the 
increase of the cross section in the ISGQR region, especially for the L=2 state whose main 
component is the ISGQR. This is at variance with the relativistic Coulomb excitation studied 
in ref. |^ because the Coulomb interaction very selectively populates dipole transitions and 
therefore cannot excite the most important two-phonon components of the ISGQR which 
are built with monopole and quadrupole phonons (see table fT|). 

The obtained ratio between cross section in the Giant Resonance region and that in the 
two phonon one varies from 3.7 in the anharmonic and non-linear case to 4.6 in the harmonic 
and linear calculation. If we only consider the cross section to the single and double isoscalar 
giant quadrupole resonance those ratios increase to 6.5 and 9.6, respectively. Those values 
are smaller than the ones reported in ref. for the cross sections at the grazing angle. This 
difference can be traced back to the present availability of the experimental elastic cross 
section needed to fix the imaginary part of the optical potential and to the fact that the 
theoretical approach has been improved in several aspects, especially in the calculation of 
the form factors. 

Our calculation can be compared with the experimental data of ref. where the re- 
action '^^Ca+^^Ca at 50 MeV/u has been studied. Let us resume the important results of 



ref. and the most critical points. We discuss first the inclusive spectrum and later on we 
will analyse the one obtained in coincidence with backward emitted particles. The inelastic 
spectrum was extracted for ejectiles scattered between 3.4 and 10 degrees in the center of 
mass frame. The GR contribution was obtained from the inclusive inelastic spectrum by 
deconvolution of the angular distributions into inelastic excitations and a non-inelastic back- 
ground. For the inelastic excitation, a DWBA prediction was used. As for the background, 
its angular distribution was assumed to be similar to the one of the energy region located 
immediately above the GRs. This procedure gave 113 mb/sr between 12 and 22 MeV for 
the inelastic excitation corresponding to 40% of the quadrupole EWSR. However, it should 
be noticed that the estimate of the non-inelastic background underlying the GR is not un- 
ambiguous. Indeed, if inelastic excitation is still present in the region above the resonance 
as expected from fig.l, the assumed background is overestimated. In this case the extracted 
value should be understood as a minimum. The maximum inelastic contribution compatible 
with the measured angular distribution is 223 mb/sr. This corresponds to the other extreme 
when no non-inelastic background is considered. Therefore the GR cross section extracted 
from the inclusive spectrum is between 113 and 223 mb/sr depending upon the background 
hypothesis. The associated EWSR would thus range between 40 and 80% if the whole cross 
section is assumed to be coming from quadrupole states. 

In order to get the total cross section one has to extrapolate the measured differential 
cross section beyond the solid angle covered by the ejectile detector. This was done by 
assuming that the DWBA angular distribution used to fit the measured angular distribution 
in ref. Q was also valid in the region in which no data are available. The ratio between the 
integrals of the DWBA cross section over the full angular range and that over the angles 
covered by the detector is 3.16. Taking into account the fraction of the solid angle covered by 
the spectrometer one gets a total compensating factor of 6.67x10"^. Such factor transforms 
the double differential cross section into the energy differential one. The resulting total 
cross section is then 7.5 and 15 mb respectively. These values have to be compared with 
the theoretical inelastic cross section which, in the anharmonic and non-linear case, adds 
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up to 22 mb in the GR region. Keeping into account the uncertainties of the analysis 
of the experimental data and the fact that our theoretical results are obtained without 
adjusting any parameter, the comparison can be considered satisfactory. In order to draw 
quantitative conclusions one should elaborate on different issues both from the experimental 
and theoretical sides. A recent experiment on the same reaction [113 using an improved 



apparatus is expected to eliminate most of the experimental uncertainties. These new data 
will allow a more reliable determination of some parameters entering in the theoretical 
calculation, mainly in the optical potential. 

Coincidences with backward emitted particles provide an unambiguous signal for the 
inelastic excitations and could in principle be used to avoid the non-inelastic background 
problems. This was the idea of ref. [Q, but some other sources of uncertainties appear. 
The coincidence rate with backward emitted protons was converted into a differential cross 
section correcting for the energy dependence of the proton multiplicity. At that time it was 
already stressed that this correction factor can be subject to many uncertainties. First of 
all, this proton multiplicity function was calculated with a statistical decay code which does 
not include any direct decay component. Furthermore, due to the absence of out-of-plane 
detectors, the azimuthal angular distribution was not measured and was assumed to be 
uniform. This procedure gives a cross section for the GR extracted from the coincidence 
data (339 mb/sr) larger than the one obtained from the inclusive inelastic spectrum (between 
113 and 223 mb/sr). This shows that the hypotheses used are not correct. The use of a 



47r detector in a recent experiment [|T^ should solve these ambiguities since it will provide 
the angular distribution of the emitted protons and there will be no need to rely on a 
statistical code to infer their multiplicity. That was not the case in the experiment of ref. [Q. 
Therefore, only the ratio was deduced from the coincidence data. Two values of this ratio 
were reported by assuming two backgrounds for the two-phonon region, while the GR peak 
was considered with no background subtraction in the coincidence spectrum. The values of 
the second phonon cross section were, after subtraction of the two backgrounds, 30 and 17 
mb/sr respectively for an energy running from 28 to 40 MeV, while the GR cross section was 
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339 mb/sr in the coincidence spectrum in the range of 12 to 22 MeV, leading to the ratios 
11 and 20 quoted in ref. Such values are the ratios between the single GR cross section 
and only a small fraction of the DGR cross section. We want to stress here that the correct 
procedure should be not to subtract any background in the two-phonon region. Indeed, 
on one hand, coincidence with backward emitted particles avoids any contribution from 
non-inelastic background in the experimental data. On the other hand, in our theoretical 
calculation, not only double GQR has been included but many contributions from different 
inelastic excitations have been taken into account. These two remarks plead in favour of a 
direct comparison of the one and the two-phonon regions with no background subtraction. 

In order to have a more direct comparison we present, in fig. 2, the experimental coin- 
cidence inelastic spectrum of ref. |^ (fig. 16 b) with no background subtraction. The right 
scale is the double differential cross section while the left scale is the energy differential cross 
section obtained with the above mentioned factor of 6.67x10"^. In the figure we present 
the theoretical results smoothed by a Lorentzian of 5 MeV width, rather than the 3 MeV 
used in fig. 1. From the figure we see that with this value the shape of the experimental 
peak in the GR region is well reproduced. It should be noticed that some contribution to 
the experimental cross section is present just below 14 MeV. However, due to the proxim- 
ity of the proton emission threshold, the correction for the multiplicity is more delicate in 
that energy region. Disregarding these two points, the overall agreement between theory 
and experiment is rather satisfactory. A rough estimate of the one-phonon and two-phonon 
cross-sections can be obtained by integrating both curves in the energy ranges shown in 
fig. 2 as shadowed areas. By doing that one would get an experimental and theoretical ratio 
of 2.4 and 2.3, respectively. We want to stress that the experimental ratio quoted above 
is different from the one deduced in ref. |Q because the latter is the ratio between the full 



^one digit inversion error was spot in the text of fig. 16 c) and d) of ref. [4] (erratum to be 
published) 
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peak of the single GR and the DGR with background subtraction, while the former one is 
obtained without background subtraction in both single GR and DGR. Furthermore the first 
two experimental points in fig. 2 were not included as explained before. Finally, we would 
like to comment on the dependence of the theoretical ratios upon the smoothing width. 
This ratio is decreasing with the increasing width, due to the fact that while the integral of 
the single GR is decreasing the one over the region of the DGR remains almost unchanged. 
This is related to the fact that in the single GR energy region the peaks of the single $q, 
states are quite separate while the density of states in the DGR region is very high. In any 
case the dependence on the width is not very strong: by varying F from 3 to 6 MeV the 
ratio changes from 2.75 to 2.20. These values cannot be directly compared with the values 
reported in tables p^[V| because the latter have been obtained just by summing the cross 
sections associated with each discrete state. 



IV. CONCLUSIONS 

We have calculated the inelastic scattering cross sections of one- and two-phonon states 
for the '^''Ca + '^°Ca collision at E/u=50 MeV. Several effects have been evidenced. In 
particular, we have analyzed the role played by anharmonicities in the excitation spectrum 
and non-linearities in the operator describing the mutual interaction of the collision partners. 
The anharmonicities are particularly important at relatively low energy where the excitation 
comes through the one-phonon component of the mixed states. The non-linearities give 
their main contribution at high energy, in particular in the region of the double quadrupole 
giant resonance. Namely, in the interval between 28 and 38 MeV, they give an increase 
of about 40% with respect to a harmonic and linear calculation. This increase is due to 
the excitation of other two-phonon states which are populated because of the presence 
of the anharmonicities and non-linearities. With all the previously discussed caveats, the 
comparison of the smoothed theoretical result with the experimental coincidence inelastic 
spectrum of ref. is satisfactory. The inclusion of three-phonon states in the calculation 
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will increase the inelastic cross section at higher excitation energies. At the same time, a 



fraction of the population of the two-phonon states will move to higher energies. In ref. |n 
it has been shown within a simple model that the spectrum calculated by diagonalizing 
in a space including up to three phonons the hamiltonian obtained by a boson expansion 
truncated at the quartic order is in reasonable agreement with the exact one. A similar 
calculation is feasible also in a realistic case. This, together with the results shown here, 
encourages us to proceed in the direction of calculating the three-phonon excitation cross 
section for the system ^''Ca + ^°Ca at E/u=50 MeV for which experiments have already 



been done [13 
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FIGURES 



FIG. 1. Inelastic cross section for the system ^°Ca + ^'^Ca at 50 MeV/u as function of the 
excitation energy. Both curves are the result of a smoothing procedure with a Lorentzian with a 
width r=3 MeV. The shadowed areas are the energy regions over which we have summed the cross 
sections reported in the tables. 



FIG. 2. The dots represent the experimental coincidence inelastic spectrum of ref. [4] (Fig. 16 
b) with no background subtraction (right scale). The solid line is the result of a smoothing 
procedure with a Lorentzian with a width r=5 MeV of the theoretical inelastic cross section 
for the anharmonic and non- linear case (left scale). The shadowed areas are the energy regions 
over which we have integrated the energy differential cross sections. The resulting values in mb 
are the numbers reported in the two areas. Those above the curves refer to the theoretical results, 
while the ones below refer to the experimental data. In the inset we report the ratios between the 
single GR cross section and the DGR ones for the two cases. 
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TABLES 

TABLE I. RPA one-phonon basis for the nucleus ^°Ca. For each state its spin, parity, isospin, 
energy and percentage of the EWSR are reported. 



Phonons 


J"" 


T 


E{MeV) 


%EWSR 


GMRi 


0+ 





18.25 


30 


GMR2 


0+ 





22.47 


54 


GDRi 


1- 


1 


17.78 


56 


GDR2 


1- 


1 


22.03 


10 


ISGQR 


2+ 





16.91 


85 


IVGQR 


2+ 


1 


29.59 


26 


3- 


3- 





4.94 


14 


LEOR 


3- 





9.71 


5 


HEOR 


3- 





31.33 


25 
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TABLE II. Characteristics of the |$q > quadrupole 2+ states whose major components are 
in the first column. In the second column we show the energies of the major components in the 
harmonic approach. The shift in the energy produced by the anharmonicities is indicated by 
A£' (in KeV). We can compare these values with the diagonal matrix elements of the residual 
interaction, A£^o (in KeV) . In the last columns we report the amplitude with which the single and 
double ISGQR components appear in the mixed states. 
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3- <i 


5 HEOR 


36.272 
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TABLE III. Coulomb plus nuclear excitation cross section for ^'^Ca + ^°Ca at 50 MeV/u. 
Each multipolarity contribution is shown for several anharmonic and non-linear combinations. 
The values for L=l and 5 are very small and they are not shown. The cross sections (in mb) are 
summed over the energy region (0 < £^ < 12 MeV). 



Phonons 


harm. & lin. 


harm. &; non-lin. 


anh. &; lin. 


anh. &; non-lin. 


L=0 


0.1 


0.3 


1.3 


2.3 


L=2 


0.2 


0.4 


0.2 


0.1 


L=3 


14.2 


16.9 


14.3 


16.8 


L=4 


0.2 


0.3 


0.2 


0.3 


L=6 


0.5 


0.7 


0.4 


0.7 


total 


15.2 


18.6 


16.4 


20.2 
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TABLE IV. Same as table |T| but for the double ISGQR region. The cross sections (in mb) are 
summed over the energy region (28 MeV < E < 38 MeV). The values in parentheses correspond 



to the double ISGQR state. 


Phonons 


harm. & lin. 


harm. &; non-lin. 


anh. &: lin. 


anh. & non-lin. 


L=0 


0.2 (0.15) 


0.3 (0.26) 


0.2 (0.14) 


0.3 (0.21) 


L=2 


0.6 (0.33) 


1.0 (0.51) 


0.6 (0.33) 


1.1 (0.53) 


L=3 


2.2 


2.5 


2.3 


2.5 


L=4 


1.0 (0.90) 


1.9 (1.83) 


0.9 (0.85) 


1.8 (1.73) 


L=6 


0.2 


0.2 


0.2 


0.2 



total 4.2 (1.38) 5.9 (2.60) 4.2 (1.32) 5.9 (2.47) 



TABLE V. Same as table [II but for the ISGQR region. The cross sections (in mb) are summed 



over the energ 


;y region (14 MeV < E < 20 MeV). In this region there are no 


states with L=3 and 


5. The values 


in parentheses 


correspond to the ISGQR state. 






Phonons 


har. & lin. 


harm. & non-lin. 


anh. & lin. 


anh. & non-lin. 


L=0 


2.7 


2.8 


2.2 


2.2 


L=l 


3.0 


2.9 


3.6 


3.3 


L=2 


13.3 (13.2) 


16.1 (16.0) 


13.8 (13.6) 


16.0 (16.0) 


L=4 


0.1 


0.1 


0.1 


0.1 


L=6 


0.2 


0.3 


0.2 


0.3 


total 


19.3 


22.2 


19.9 


21.9 
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E (MeV) 




E (MeV) 



